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Abstract 

3 P 2 pairing in neutron matter is investigated using the Bonn potential models. We 
find pairing energy gaps in pure neutron matter comparable to the results of previous 
investigators when the attractive tensor coupling is included. However, taking into 
account that in a neutron star we have matter at j3 equilibrium, we find that the 
3 P 2 - 3 F 2 energy gap is reduced considerably. 



1 Introduction 



Superfluidity in neutron stars is expected to have a profound influence on their 
dynamical and thermal properties [1]. Due to the properties of nuclear forces 
several nucleon superfluids can be formed, the most commonly considered 
being the l So neutron and proton superfluids, and the 3 P 2 neutron superfluid. 
In a recent investigation [2] we calculated 1 S'o neutron and proton pairing 
gaps, ignoring complicated many-body effects. In this paper we follow the 
same philosophy in calculating energy gaps for neutrons in the 3 P 2 state, 
using realistic nucleon-nucleon (NN) forces derived from meson theory [3]. 



2 Theory 



The standard BCS theory [4,5] has to be extended when dealing with pairing 
of particles in states with non-zero angular momentum. The method of gen- 
eralizing the Bogolyubov transformation to pairing of particles with arbitrary 
angular momentum and isospin is described in refs. [6,7]. We are in this paper 
concerned with neutron pairing in the 3 P 2 - 3 F 2 state, and thus take J = 2, 
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S = 1, L = 1, 3, and T = 1. In this case, the l So gap equation is replaced by a 
set of ten (2 x (2J+ 1)) coupled equations for the ten energy gap components: 
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where A = (T = 1, 5" = 1, J = 2), the orbital angular momenta L, L', and L" 
take on the values 1 and 3, and 

G ALmj (k) = { ^<7i^<72 1 Sms^ (Sm s Lm L \Jmj)Y Lmj - ms (k)^ 



is a 2x2 matrix in spin-^ space. Furthermore, Yhmj-m s is a spherical harmonic, 
D\ is given by 

^(k) = ^E E A ALmj (A;)A ALm:/ (A;)Tr[GL mj (k)G AL , m;/ (k)], 

LL' mjm'j 



and efc and are single-particles, evaluated at k = k and k = kp, respectively. 
The symbol Tr denotes the trace operation, the symbol ' denotes Hermitian 
conjugation, and (k'\V^ L ' \k^j are matrix elements of the NN interaction in 
momentum space, in the partial wave characterized by L, L' and A. 



2.1 Uncoupled equations for 3 P 2 pairing 



The set of equations (1) is rarely solved in its full generality, and we will in this 
paper look at two approximations. In the first one, considered by e.g. Baldo 
et al. [8], the tensor interaction is neglected, and thus pairing in a pure 3 P 2 
wave is considered. Furthermore, only the solutions with mj = or \mj\ = 2 
are considered, because these have been found to be close to the most general 
solution of eq. (1) [7]. It can then be shown that the gap equations reduce to 

Mk) = ~- ]dk>k»(k>\V( 3 P 2 )\k)£M, < = 0,2, (2) 
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where % — 0, 2 for rrtj = 0, |mj| = 2 respectively. We have here defined 
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where -E'(k) is the quasiparticle energy, given by 
£i(k) 



e fc -e, F ) 2 + ^A 2 (l-cos 2 £), i = 2, 
^( efe - efcF )2 + _l_A2(fc)(l + 3cos 2 0), 2 = 0. 



(4) 



In eq. (2) (/c / |V r ( 3 P 2 )|&) are the momentum-space matrix elements of the bare 
NN interaction in the 3 P 2 wave. Eq. (2) must in each case be solved self- 
consistently for the gap component A (k) or A 2 (k). The energy gap is by 
convention defined to be the second term appearing under the square root in 
the expression for the quasiparticle energy: 
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When giving numerical results, we use angle-averaged energy gaps, given by 



Dl{k) = -Aj(k) 

Replacing angle- dependent quantities in the gap equation with their angular 
average has been found to be a good approximation [8]. However, the angle 
dependence is easily handled numerically, so we have throughout solved the 
gap equation with angle dependencies included. 



2.2 Coupled equations for 3 P 2 - 3 F 2 pairing 



In the second approximation we include the tensor coupling effect, but look 
at the solution with mj = only, thus reducing the number of components 
from ten to two. This solution is expected to be a good approximation to the 
most general solution of eq. (1) [7]. For mj = 0, the two coupled equations for 
the L = 1 and L = 3 gap components become 
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The functions of the polar angle of k' are defined by 
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the quasiparticle energy E(k) is given by 

E(k) = ^(e k -e kF r + Dl(k), (9) 

and the energy gap by 

D 2 x (k) = i/(0)A 2 lo (A:) + ^(0)A 2 3O (AO + g(9)A xw (k)A X30 (k). (10) 

In ref. [8] Baldo et al. obtained a large 3 P 2 energy gap without the tensor 
coupling. However, in refs. [7,9] the added attraction from the tensor coupling 
effect was found to be essential to the existence of superfluidity in this state. 
The importance of the tensor coupling depends to some extent on the choice 
of NN interaction, and we will therefore consider both cases. 

2. 3 Calculation of single-particle energies 



The single-particle energies appearing in the quasiparticle energies (4) and 
(9) are obtained through a self-consistent BHF calculation, using a G-matrix 
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defined through the Bethe-Brueckner-Goldstone equation as 

G = V + V Q TT G, (11) 



where V is the nucleon-nucleon potential, Q is the Pauli operator which pre- 
vents scattering into intermediate states prohibited by the Pauli principle, H 
is the unperturbed Hamiltonian acting on the intermediate states and u is the 
so-called starting energy, the unperturbed energy of the interacting particles. 
Methods to solve this equation are reviewed in ref. [10]. The single-particle 
energy for state ki (i encompasses all relevant quantum numbers like momen- 
tum, isospin projection, spin etc.) in nuclear matter is assumed to have the 
simple quadratic formQ 



where m* is the effective mass. The terms m* and 5, the latter being an 
effective single-particle potential related to the G-matrix, are obtained through 
the self-consistent BHF procedure. The so-called model-space BHF method 
for the single-particle spectrum has been used, see e.g. ref. [10], with a cutoff 
momentum ku = 3.0 fm -1 > k F . In this approach the single-particle spectrum 
is defined by 

k 2 

e k . = — + Ui, (13) 



where m is the nucleon mass, and the single-particle potential Ui is given by 

, E kh <k F (kik h \G(uj = e ki + e kh )\kik h ) , k { < k M , 
Ui = \ , (14) 

0, h > k M , 



where the subscript AS denotes antisymmetrized matrix elements. This pre- 
scription reduces the discontinuity in the single-particle spectrum as compared 
with the standard BHF choice ku = kp. The self-consistency scheme consists 
in choosing adequate initial values of the effective mass and S. The obtained 
G-matrix is in then used to calculate the single-particle potential from 
which we obtain new values for m* and 5. This procedure continues until 
these parameters vary little. We stress that in the energy gap calculations we 
have used the single-particle spectrum defined by eqs. (13) and (14), and not 
the effective mass approximation (12). 

1 Throughout this work we set H = c = 1. 
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The BHF equations can be solved for different proton fractions, using the 
formalism of refs. [10,11]. The conditions for f3 equilibrium require that 

/J, n = Hp + He, (15) 



where Hi is the chemical potential of particle type i, and that charge is con- 
served 

n p = n e , (16) 



where rii is the particle number density for particle i. If muons are present, 
the condition for charge conservation becomes 

n p = n e + n fl , (17) 



and conservation of energy requires that 

He = /V ( 18 ) 



We assume that neutrinos escape freely from the neutron star. The proton 
and neutron chemical potentials are determined from the energy per baryon, 
calculated self-consistently in the MBHF approach. The electron chemical po- 
tential, and thereby the muon chemical potential, is then given by He = Hn—Hv 
The Fermi momentum of lepton type I — e, h is found from 

k Fl = h] -rnf (19) 



where mi is the mass of lepton /, and we get the particle density using 
ni = kf/37r 2 . The proton fraction is then determined by the charge neutrality 
condition (17). 

Throughout this work, we use the meson-exchange potential models of the 
Bonn group, versions A, B and C in Table A.l of ref. [3]. These potentials 
have given satisfactory results in applications to finite nuclei and nuclear mat- 
ter [10]. Furthermore, the only parameters entering in the Bonn models are 
physically motivated ones, e.g. meson masses. A point of interest in this work 
is that the Bonn potentials have a rather weak tensor component. The A.l 
parameters are the consistent choice of parameters in a non-relativistic ap- 
proach, as they were fitted using the Blankenbecler-Sugar reduction [12] of 
the four- dimensional Bethe-Salpeter equation, yielding a T-matrix equation 
of the same form as the non-relativistic Lippmann-Schwinger equation [3]. 
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3 Pairing in the pure 3 P 2 wave 



We have solved the gap equations (2) in pure neutron matter. The calculational 
procedure consists of two steps: First we obtain self-consistent single-particle 
energies in the MBHF approach, as described in section 2.3. The next step 
is to solve the relevant gap equation, for mj = or \mj\ = 2, with the bare 
NN interaction as pairing interaction. Being non-linear integral equations, the 
gap equations must be solved by iteration. In ref. [2] we applied a model- 
space approach to the gap equation for 1 S'o pairing. One of the advantages of 
this approach was that the numerical treatment became easier, as it allowed 
us to treat pairing correlations and the effects of the repulsive core of the 
NN interaction separately. The same method can easily be applied to the gap 
equations studied in this paper. However, as seen from fig. 1, the interaction in 
the 3 P 2 wave is weak at the densities of interest, so a model space treatment 
is not called for. Nevertheless, we would like to point out some problems 
encountered in the numerical solution. First, the 3 P 2 energy gap is considerably 
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Fig. 1. Matrix elements (kp\V( 3 P2)\k) of the Bonn A, B and C potentials, shown 
as functions of the relative momentum k for kp = 1.8 fm -1 . 

smaller than the 1 So energy gap, which is of order 1 MeV [2]. A small energy 
gap at the Fermi surface causes problems in the numerical treatment, as we 
see from eq. (2) that for small values of Aj(fcp) the integrand 

mk% = -h-(k\vep 2 )\k')^M , = o,2 
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in the gap equation has a very narrow peak at k — /cp, as shown in fig. 2. 
Thus, to obtain stab le results, we had to take great care in choosing mesh 




k (fm ) 



Fig. 2. Integrand f(k, k')o in the gap equation for mj = and k = kp = 2.4 fm 



points for the momentum space integration near kp. As seen from fig. 2 the 
choice of upper cutoff momentum was almost immaterial, and k — 15 fm -1 was 
found to be adequate. We should mention that Khodel, Khodel and Clark [13] 
recently published a new method for solving gap equations which eliminates 
the problems occuring when the gap parameter is small. 

In fig. 3 the energy gap D (kp), calculated with MBHF single-particle energies, 
is shown. The energy gap Z^&f) was found to be similar to Do(kp), as seen 
from fig. 4. Takatsuka and Tamagaki [7] and Amundsen and 0stgaard [9] 
found that the 3 P 2 gap vanished for realistic values of the neutron effective 
mass (in the range 0.90-0.65, see table 1), and thus it was necessary to include 
the attractive tensor coupling to the 3 F 2 wave. As will be shown in the next 
section, our results for the pure 3 P 2 wave are of the same size as the results 
where the tensor coupling is included. This result is probably due to the Bonn 
potentials having a weaker tensor force than the potentials used in refs. [7,9]. 
However, our results are one order of magnitude smaller than the results of 
Baldo et al. [8] , who found energy gaps of the order 1 MeV without the tensor 
coupling effect. 

The energy gap function D given by eq. (5) is nodeless (as a function of 
0), while D 2 has nodes for cos 2 = 1. A nodeless energy gap results in the 
usual exponential decay of the specific heat below the superfluid transition 
temperature. If there are nodes, the specific heat will follow a power-law [14]. 
Thus the solution with m,j = 0, if realized, gives a different thermal behaviour 
than the solution with \mj\ = 2. The favoured solution is the one which gives 
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Fig. 3. Angle-averaged energy gap Do(kp) calculated with Bonn A, B and C. 



o : m = 




k p (fm 1 ) 

Fig. 4. Comparison of the angle averaged energy gap for mj = with corresponding 
results for \mj\ = 2. The results are for the Bonn B potential. 

the lowest condensation energy. This quantity is, in the so-called weak-coupling 
limit, given by 

AE « - P j-D\{k F \ (20) 

where po is the density of states at the Fermi energy. In table 1 the two 
solutions are compared with respect to the condensation energy. The results 
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k F (far 1 ) -AE/(p /87r){mj = 0) (MeV) -AE/(po/8ir)(\mj\ = 2) (MeV) 



i An 
1.4U 


U.UUo 


n one 
U.UUo 


1.60 


0.058 


0.056 


1.80 


0.113 


0.116 


2.00 


0.139 


0.133 


2.20 


0.100 


0.095 



Table 1 



Condensation energies calculated with the Bonn B potential. 

shown are for the Bonn B potential. It is seen that the two solutions are very 
nearly degenerate, a result which is in agreement with the more complete 
discussions found in refs. [7,9], with the nodeless solution, mj = 0, seemingly 
the favoured one. 



4 Inclusion of the tensor coupling 

We have solved the coupled equations (7) and (8) for pure neutron matter. 
We give results for the angle average of D\, 




where the last approximation is valid because Aa3o(&f) was found to be negli- 
gible in comparison with Aaio(^f)- The following procedure was used to solve 
the coupled gap equations: First we obtained single-particle energies in the 
MBHF approach, as described in section 2.3. Starting with suitable approxi- 
mations for the gap components A A10 and A A30 (for A A10 we could start with 
the solution of the uncoupled equation considered in section 3), the gap equa- 
tions (7) and (8) were solved by iteration. Again, special attention had to be 
paid to the integration over momenta near fop, because of the almost singular 
behaviour of the integrand there. Stable results were obtained with the same 
choice of mesh points as for the uncoupled case considered in the previous 
section. 

In fig. 5 the energy gaps in pure neutron matter obtained with Bonn A, B 
and C are shown. From these results it is seen that the tensor force increases 
the energy gaps. However, this effect is not so important here as in the works 
of Takatsuka and Tamagaki [7] and Amundsen and 0stgaard [9]. This was to 
be expected, since the Bonn potentials have weaker tensor components than 
the OPEG potentials used in refs. [7,9]. The effects of the tensor coupling 
are greater for Bonn C than for Bonn B and A, reflecting the differing tensor 
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Fig. 5. 3 P2- 3 -?2 energy gaps in pure neutron matter. 

strengths of the Bonn interactions. In fig. 6 the effective masses are plotted as 
functions of the neutron Fermi momentum. This figure shows that the Bonn 
models give nearly identical effective masses in pure neutron matter. The 
differences in the energy gaps are therefore mainly due to differences in the 
3 P 2 - 3 F 2 wave of the interactions. From fig. 1 it is seen that Bonn C gives more 
attraction near kp than do Bonn B and A. Since the main contribution to the 
integral in the gap equation comes from the region near kp, it is reasonable 
that Bonn C gives the largest energy gaps. 
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Fig. 6. Effective masses in pure neutron matter and in matter at (3 equilibrium. 
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We have also calculated the 3 P2- 3 -^2 neutron energy gap in neutron star matter 
at (3 equilibrium. The neutron fractions, shown in figure 7, were determined 
from a self-consistent MBHF calculation, imposing the relevant equilibrium 
constraints as described in section 2.3. The results, obtained with the Bonn 
A potential, are shown in figure 8. The size of the energy gap is reduced 



0.00 0.10 0.20 0.30 0.40 0.50 

n (fm~ 3 ) 

Fig. 7. Neutron fraction as a function of the nucleon density n at (3 equilibrium with 
muons included. 



LU 0.040 - 




k p (fm ) 



Fig. 8. Neutron energy gap at (5 equilibrium. 

considerably, and the density range for this superfluid shrinks in comparison 
with the case of pure neutron matter. The main reason for this is probably 
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the reduced neutron effective masses in /3-stable matter, as seen in fig. 6. 



5 Discussion 



The results of this paper are summarized in fig. 9, where we show the following 
results: 

• the 3 P2- 3 -^2 energy gap in pure neutron matter 

• neutron 3 P2- 3 P2 energy gap at f3 equilibrium 

• proton 1 5'o energy gap at j3 equilibrium, taken from ref. [2], 

all as functions of the nucleon density. From the figure it is seen that the 
proton energy gap is one order of magnitude larger than the neutron energy 
gaps. Thus the dominating contribution to the suppression of the modified 
URCA processes will be due to superconducting protons [15]. Taking into 
account that in a neutron star we have matter at j3 equilibrium and not pure 
neutron matter, the neutron energy gap is reduced by approximately a factor 
of two. 



> 0.6 - 




0.00 0.10 0.20 0.30 0.40 0.50 

n (fm ) 

Fig. 9. Energy gaps for protons and neutrons in the quantum fluid region of a 
neutron star. 



The results obtained in this work for the 3 P2 _3 ^2 energy gaps in pure neutron 
matter are in good agreement with those of Takatsuka and Tamagaki [7] and 
Amundsen and 0stgaard [9]. We find a maximum energy gap of the same size 
as in refs. [7,9], and a density region for pairing intermediate between those 
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Authors 


Max. energy gap (MeV) 


Density range (fm 1 ) 


m*/m (typical values) 


Potential 


A0 


0.12 


1.4-3.0 


0.70-0.89 


OPEG 


BCLL 


1.5 


1.8-3.2 


0.70-0.86 


Argonne V\4 


TT 


0.12 


1.5-2.4 


0.70-0.82 


OPEG 


HGRR 


0.9 


1.0-3.0 


1.0 


Tabakin 


This work 


0.13 


1.2-2.5 


0.65-0.91 


Bonn B 



Table 2 

Results for the 3 P 2 energy gap in neutron matter. 



found in the two references above. In table 2 and fig. 10 we summarize the re- 
sults of refs. [7](TT), [8] (BCLL), [9](A0), [16] (HGRR), and the present work. 
We should at this^oint emphasize that the effective masses quoted in the 

o : Bonn B 

[]: Amundsen & Oestgaard 
+ : Baldo et al. 




1.0 1.5 2.0 2.5 3.0 

k F (fm" 1 ) 



Fig. 10. Results for the triplet energy gap in neutron matter. The results of Baldo 
et al. and Hoffberg et al. were obtained without the tensor coupling. 

table are meant to facilitate comparison between the different works. Both in 
this work and in refs. [7,8], the gap equations were solved with fc-dependent 
single-particle potentials, without resorting to the effective mass approxima- 
tion. From table 2 it is seen that the the energy gap has generally been found 
to be of order 0.1 MeV. The results of Baldo et al. [8] are remarkable, con- 
sidering that they have neglected the attractive contributions from the tensor 
coupling to 3 F 2 . With effective masses in the range used by Baldo et. al, the 
energy gap in the pure 3 P 2 wave was in this work found to be one order of 
magnitude smaller, and in refs. [7,9] it was found to vanish. In addition, we 
found the 3 P 2 gap to exist at lower densities than Baldo et al. Although both 
Amundsen and 0stgaard and Takatsuka and Tamagaki used the OPEG po- 
tentials, Amundsen and 0stgaard found that the energy gap extended up to 
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kp — 3.0 fm -1 . This difference can be understood because different approaches 
to the single-particle energies were employed. Amundsen and 0stgaard took 
m*/m to be constant (equal to 0.70) above kp = 2.4 fm -1 [18], while Takat- 
suka and Tamagaki used a density-dependent effective mass all the way, and 
thus it dropped to values below 0.70 beyond kp = 2.4 fm -1 . Thus it is seen 
that the single-particle spectrum is a crucial ingredient in energy gap calcula- 
tions. However, the Brueckner-Hartree-Fock spectrum, and variations thereof, 
used in most works, including the present, may be of questionable validity at 
the densities where the 3 P2 energy gap closes, where relativistic effects and 
many-body correlations become increasingly important. To further illustrate 
the importance of the single-particle energies, we have calculated the energy 
gap at kp = 1.8 fm -1 as a function of the effective mass parameter m*/m. The 
results in fig. 11 show that the energy gap depends strongly on the effective 
mass. Also include^dj n table 2 are the results of Hoffberg et a l. [16]. They are, 



0.70 - 




0.00 1 ' 1 ' 1 ' 1 ' 1 

0.70 0.80 0.90 1.00 1.10 

m*/m 

Fig. 11. Effective mass dependency of the 3 P2- 3 P2 energy gap. The results were 
obtained using the Bonn C potential and kp = 1.8 fm -1 . 

however, not directly comparable with the other results, as they have been 
obtained directly from the 3 P2 phase shifts, and with m*/m = 1. However, it 
is interesting to note that they obtained a smaller maximum energy gap than 
Baldo et al. 

To the best of our knowledge, the 3 P 2 - 3 P 2 energy gap has not been calculated 
in /3-stable neutron star matter before. Our results show that the energy gap 
is reduced by a factor of two when (3 equilibrium is taken into account, mainly 
because the neutron effective mass is reduced. A point of interest here is 
that in a recent investigation of neutron star cooling, Schaab et al. [19] found 
that agreement with observations was obtained if the 3 P 2 energy gaps were 
lowered by a factor of two (they used the results of ref. [9] for pure neutron 
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matter). A recent analysis by Page [20] of data for the Geminga pulsar showed 
that baryon pairing in most of the star was necessary to obtain agreement 
with the observed surface temperature, assuming that fast cooling processes 
were possible. However, with equations of state for /3-stable matter, calculated 
within the non-relativistic BHF approach, we find that the critical density 
for the onset of the direct URCA processes to be 0.62 fm -3 , while within the 
relativistic Dirac-Brueckner-Hartree-Fock approach, the corresponding density 
was found to be 0.50 fm~ 3 [15]. At these densities the neutron and proton 
pairing gaps have already vanished, so the direct URCA processes will not 
be suppressed by superfluidity within this picture. Further investigation of 
3 P2 pairing is therefore of interest. However, an important shortcoming of all 
calculations, including those reported in this work, is the neglect of many- 
body effects in the pairing interaction. The so-called polarization term has 
been found to reduce the 1 Sq proton and neutron energy gaps by as much as a 
factor of 2-3 [21]. It is, however, expected that polarization effects will tend to 
increase the 3 P2 neutron energy gap [22]. Further investigation of superfluidity 
in dense matter is therefore of interest. 



6 Conclusion 

In this work we have calculated pairing gaps for neutrons in the 3 P2 state in 
neutron matter and in matter at (3 equilibrium. For pure neutron matter, we 
found a non-zero energy gap for 1.2 fm -1 < fop < 2.5 fm -1 , with a maximum 
value of 0.12 — 0.13 MeV. Taking the conditions for (3 equilibrium into account, 
the corresponding values became 1.32 fm -1 <kp< 2.34 fm -1 and 0.05 MeV. 
However, evaluation of the contribution from medium polarization effects to 
the pairing interaction is necessary before any definite conclusions can be 
drawn. 

This work has been supported by the Research Council of Norway (NFR) 
under the Programme for Super computing through a grant of computing time. 
MHJ thanks the NFR and the Istituto Trentino di Cultura, Italy, for financial 
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